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The  formulation  of  expressions  for  transport  prop¬ 
erties  in  nonequilibrium  flows  has  been  an  active 
research  field  for  many  years  [?].  One  powerful 
feature  of  particle  simulation  methods  such  as  di¬ 
rect  simulation  Monte  Carlo  (DSMC)  is  that  they 
do  not  require  transport  properties  as  input  param¬ 
eters.  Rather,  given  a  sufficiently  realistic  model 
of  the  intermolecular  potential  and  energy  trans¬ 
fer,  transport  properties  emerge  naturally  as  a  sta¬ 
tistical  consequence  of  many  collisions  along  with 
boundary  conditions  [?].  Thus,  a  physically  realis¬ 
tic  yet  computationally  tractable  model  for  molec¬ 
ular  collisions  is  of  primary  importance  for  DSMC 
methods. 

In  a  DSMC  simulation,  collision  pairs  are  se¬ 
lected  based  on  the  local  density  and  the  velocity- 
dependent  collision  cross  section.  For  any  given 
reaction  (or  relaxation)  model,  the  number  of  re¬ 
actions  produced  in  a  given  calculation  (the  local 
reaction  rate)  will  be  controlled  by  the  local  to¬ 
tal  collision  frequency.  Therefore,  any  increase  or 
decrease  in  a  molecules  collision  cross  section  can 
directly  affect  its  effective  reaction  rate. 

Gorbachev  et  al.  [?]  have  presented  a  derivation  of 
analytical  expressions  for  the  average  internuclear 
distance  R(v.J)  of  a  diatomic  molecule  as  a  func¬ 
tion  of  vibration  and  rotation  level,  based  on  an  ac¬ 
curate  internuclear  potential.  Gimelshein  et  al.  [?] 
have  presented  an  implementation  of  these  expres¬ 
sions  in  a  DSMC  code  as  an  addition  to  the  variable 
soft  sphere  (VSS)  collision  model.  This  intriguing 
proposal  will  produce  significantly  increased  colli¬ 
sion  cross  sections  for  particles  with  high  levels  of 
internal  excitation,  thus  increasing  their  rates  of  re¬ 
action.  Although  the  expressions  for  average  inter¬ 
nuclear  distance  R(v,J)  may  be  expected  to  be  quite 
accurate,  the  effect  on  intermolecular  collision  cross 
section  may  not  be  straightforward.  The  goal  of 
the  ^present  discussion  is  to  examine  the  feasibility 
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of  validating  the  realism  of  this  approach. 

Transport  properties  are  an  important  source  of  in¬ 
formation  on  intermolecular  potentials.  The  follow¬ 
ing  discussion  refers  to  viscosity,  but  most  of  the 
same  arguments  apply  to  diffusion.  The  viscosity 
is  typically  dominated  by  elastic  collisions,  so  that, 
in  a  first  approximation,  molecules  may  be  viewed 
as  atom-like.  If  we  consider  higher-order  effects, 
we  must  include  inelastic  collisions.  The  effect  of 
inelastic  collisions  on  the  transport  collision  inte¬ 
gral  (for  viscosity)  is  expected  to  be  small  (Mason- 
Monchick  approximation),  but  has  not  been  thor¬ 
oughly  investigated  for  high  temperatures.  The 
other  effect  specific  to  molecules  is  the  subject  of 
this  paper:  that  is.  higher  rovibrational  (v.J)  states 
will  increase  R(v,  J),  which  may  in  turn  increase  the 
collision  cross  section.  Since  the  population  of  high 
(v,J)  states  will  become  significant  at  high  temper¬ 
atures.  one  might  expect  that  the  viscosity  at  high 
temperatures  may  reflect  this  effect  [?].  However, 
in  addition  to  the  difficulty  in  obtaining  accurate 
viscosity  data  at  very  high  temperatures,  any  ex¬ 
amination  of  these  data  to  glean  insight  into  the 
effect  of  R(v,J)  would  need  also  to  disentangle  the 
effect  of  inelastic  collisions. 

To  estimate  the  feasibility  of  validating  the  effect 
of  increasing  R(v.J)  at  higher  temperatures  from 
viscosity  data,  some  estimates  are  provided  for  a 
simple  gas  of  pure  molecular  hydrogen.  The  (v.J) 
populations  are  in  equilibrium  and  the  viscosity  col¬ 
lision  cross  section  is  given  by  the  VSS  cross  section 
where  each  (v.J)  state  has  a  different  reference  di¬ 
ameter  as  defined  in  [4].  Fig.  1  shows  that  the  as¬ 
sumed  effect  of  R(v,J)  will  begin  to  significantly  de¬ 
crease  the  viscosity  compared  with  the  VSS  model 
for  very  high  temperatures.  However,  the  effects 
of  dissociation  and  ionization  at  these  high  tem¬ 
peratures  in  a  real  gas  are  expected  to  be  of  greater 
importance  to  the  measurable  viscosity  than  the  di¬ 
ameter  effect.  The  calculation  indicates  that  the  ef¬ 
fect  of  (v,J)  excitation  on  collision  cross  section  due 
solely  to  the  diameter  as  proposed  in  [4]  is  likely  im- 


possible  to  verify  through  viscosity  data.  However, 
other  effects  such  as  a  change  in  the  attractive  well 
depth  or  inelastic  collisions  may  contribute  in  real¬ 
ity. 


Figure  1:  Variation  in  viscosity  due  to  addition  of 
rovibrational  diameters  to  VSS  model 

The  viscosity  effect  shown  above  may  seem  to  in¬ 
dicate  that  including  collision  diameters  based  on 
(v.J)  level  is  of  no  importance.  This  conclusion 
is  not  necessarily  true  for  highly  nonequilibrium 
gases,  which  are  naturally  the  regime  of  interest  for 
DSMC  simulations.  Certain  nonequilibrium  gases 
have  highly  excited  vibrational  distributions  and,  in 
these  cases,  an  increase  in  collision  diameter  for  the 
molecules  with  high  vibrational  quantum  number 
could  change  the  chemistry  prediction  for  a  DSMC 
simulation. 

A  number  of  experiments  using  light-induced  drift 
(LID)  have  been  performed  which  can  directly  mea¬ 
sure  the  effect  of  changing  (v.J)  on  the  diffusion 
cross  section  of  certain  molecules  [?].  These 
measurements  show  that  the  diffusion  collision  fre¬ 
quency,  and  thus  the  collision  cross  section,  typi¬ 
cally  increases  a  small  amount  when  v  is  increased 
from  0  to  1.  On  the  nthpr  hand  these  molecules 
-showrrdecrease  in  collision  frequency  byaffey  per¬ 
cent  as  J  increases.  It  is  not  clear  how  similar  the 
results  would  be  for  other  molecules.  It  is  also  dif¬ 
ficult  to  assess  how  the  results  for  a  small  range  of 
(v,J)  values  would  change  for  very  high  v  or  J  lev¬ 
els.  These  results  do  demonstrate,  however,  that 
even  though  R(v,J )  increases  in  a  predictable  way 


with  increasing  internal  energy,  the  effect  on  the 
intermolecular  collision  cross  section  is  difficult  to 
predict  and  is  influenced  by  subtle  details  of  the 
intermolecular  potential. 

The  proposal  to  include  increased  collision  cross  sec¬ 
tions  for  high  (v,J)  states  in  DSMC  simulations 
is  worthy  of  examination.  It  is  based  on  a  first- 
principles  approach  to  fundamental  molecular  prop¬ 
erties.  While  the  majority  of  flowfields,  where  the 
populations  of  very  high  internal  energy  states  are 
insignificant,  need  not  consider  such  an  effect,  there 
are  certain  nonequilibrium  cases  where  this  effect 
could  potentially  be  important.  It  seems,  however, 
that  the  collision  cross  section  may  not  always  in¬ 
crease  with  (v.J)  in  the  straightforward  manner  pro¬ 
posed.  Details  of  specific  characteristics  of  the  ap¬ 
plicable  intermolecular  potential  may  need  to  be 
considered.  Some  validation  of  the  relationship 
between  internuclear  diameter  and  intermolecular 
collision  cross  section  is  recommended  before  it  is 
widely  applied  in  DSMC  codes. 
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